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SUMMARY 

The ideas of Benney & Bergeron (1969) and Davies (1970) on nonlinear critical lay- 
ers are extended, and some new nonlinear neutral modes are computed for compressible 
boundary layer flow. A special case of the work is when the generalized inflexion point 
criterion holds. Neutral modes are found for a range of phase-speeds, dependent on the 
Mach number, and the properties of these are discussed. As in the linear case when the 
flow is relatively supersonic, multiple neutral modes exist. The behaviour of the neutral 
amplitude in some limiting cases is also considered, and it is found that the results are 
significantly different from that in incompressible flow when the flow is locally 
supersonic. 
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1. Introduction 

The resurgence of interest in high speed flows has stimulated further work in the area 
of compressible boundary layer flows. Most of the early work on the linear instability 
of compressible boundary layers was carried out by Lees & Lin (1946), Lees & Reshot ko 
(1960), Reshot ko (I960) and this work and subsequent contributions are described in the 
reviews by Reshotko (1976), Mack (1984,1986). More recent work based on more formal 
asymptotic techniques includes that of Smith (1989), Gajjar & Cole (1989a), Cowley &z 
Hall (1988) for boundary layers, Goldstein & Leib (1989) for shear layers, as well as several 
other related papers dealing with various aspects of the numerical solution of the stability 
equations, see for example Malik (1987), and Jackson & Grosch (1989). It is difficult to 
conduct experiments on high speed flows and to date there are very few reliable results 
and experimental data for transition in supersonic and hypersonic boundary layers. Not 
much work has been done to address the nonlinear instability of compressible boundary 
layers, theoretical, experimental or even numerical. 

The current work forms a continuation of some earlier work described in Gajjar & 
Cole (1989a) but with the emphasis now on studying the effect of a strongly nonlinear 
critical layer on the stability on a compressible boundary layer. The initial ideas on 
strongly nonlinear critical layers were first put forward by Benney &; Bergeron (1969) and 
Davies (1970), hereafter referred to as BBD, and were subsequently applied and extended 
by several authors. An excellent account of much of this work is given by Stewartson 
(1981). The underlying key ideas may be summarized as follows. If one considers the 
linear inviscid instability of an incompressible parallel fluid flow, say Ub(Y)j to normal 
mode disturbances proportional to Se lQ ^ x ~ ct \ then the governing equation for instability 
is Rayleigh's equation. Here Y is the normal coordinate, x the streamwise coordinate, 
i the time, a,c, the wavenumber and phase-speed of the disturbance wave respectively, 
and S is the small disturbance amplitude. If the base flow is taken to be Plane Poiseuille, 
or Blasius flow with the parallel flow approximation, then there are no unstable linear 
inviscid modes. For neutral modes, since Rayleigh’s equation has a regular singularity at 
the critical point Y = Y c where Ub{Yc) = the jump in the Reynolds stress across the 
critical point has to be zero. This jump however is also proportional to the jump in the 
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logarithmic term arising in the local Frobenius expansion, across the critical point. The 
structure of the upper-branch neutral Tollmein-Schlichting mode for the flows mentioned 
above is equivalent to considering long wave limits of the Rayleigh equation with the 
inclusion of small viscous effects in appropriate thin regions. For these linear modes this 
jump is nonzero. However BBD postulated some new nonlinear neutral modes in which 
this jump was zero, and gave a description of the structure of this strongly nonlinear 
critical layer. The case of weaker nonlinearity was later considered by Haberman (1972) 
and he showed how the gap between the linear long wave upper-branch neutral modes 
and the BBD modes, could be bridged by setting up a structure in which the jump in 
the logarithmic term varied from its linear nonzero value, to zero for the BBD case. This 
jump depends on the size of the disturbance amplitude 6. The various shortcomings and 
modifications of the early work are described in Stewartson (1981). Also many of the key 
ideas were later repeated/reformulated by Smith & Bodonyi (1982), and Bodonyi et al. 
(1983), and they showed how this work could be put into a more formal asymptotic context 
as far as the stability of the flow was concerned. 

For compressible boundary layer flows Gajjar & Cole (1989a) (and hereafter referred 
to as GC), showed how the work of Haberman (1972), and Smith & Bodonyi (1982) could 
be extended to include the effects of compressibility. One of the limits considered in this 
paper was the case when the jump in the logarithmic term approached zero corresponding 
to the strongly nonlinear critical layer. The present paper considers this latter problem 
and extends the work of Bodonyi, Smith & Gajjar (1983), (and hereafter referred to as 
BSG), to compressible flows. 

The role/importance of the nonlinear critical layer with respect to the stability of the 
flow', has been considered by the authors already cited. However a lot of this work has been 
based on the assumptions of neutral or equilibrium states only. The conclusions arising 
from these studies are thus limited and can only be interpreted as giving a qualitative 
description of the downstream development of a disturbance travelling wave. These studies 
suggest a shortening of scales and the outward movement of the critical layer from the wall 
into the free stream with increased disturbance amplitudes. The important issue of w hether 
such states are stable or not cannot be deduced from this type of approach. The papers 
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by Gajjar k Smith (1985), Goldstein et ai (1987), for boundary layers, and Goldstein 
k Hultgren (1988), Goldstein k Leib (1989) together with many other similar papers, for 
shear layers, have attempted to address this very issue. From these studies it is apparent 
that the neutral states considered previously are unstable in some cases, although for shear 
layers saturated states do exist for some parameter ranges. Also despite many similarities 
in the disturbance structures and solutions, there are one or two subtle differences which 
can be important, see Goldstein k Hultgren (1988). 

In the present paper we concern ourselves with neutral states only. How relevant the 
work is is an open question since the conclusions from studies of incompressible flow may 
not necessarily carry over to those for compressible flow. In addition it is anticipated that 
the ideas and structure of the nonlinear compressible critical layer studied here will hold 
for the more general problem. Work on the full problem of following the spatial/temporal 
development of an unstable mode is in progress. 

In section 2 below the basic problem for the strongly nonlinear critical layer is formu- 
lated. Section 3 describes the results and conclusions. Much of the detailed analysis and 
implications of of the critical layer structure and the solutions is omitted here. A fuller 
account of the work may be found in Gajjar k Cole (1989b). 



2. Problem Formulation 

We adopt the same notation and nondimensionalization used in GC so that (x y y y t) 
are the nondimensional cartesian coordinates and time, (u,r) the corresponding velocities, 
(p,p, T) the pressure, density and temperat ure respectively. The local Mach number is 
A/oo , and Re is the Reynolds number which is assumed to be large. We also introduce the 
small parameter t - Re~ 1 / 12 . The basic boundary layer flow is given by, 

{u,v,p,p,T) = (Ub,i 6 Vb,Pb,pb,T b ) + 

and we assume that Chapman’s viscosity law holds. 

The disturbance structure for the flow is as sketched in Fig. 1 see also BSG, with 
I denoting the mainly inviscid region of thickness 0(e 6 ), and Y = e~ 6 y, is the boundary 
layer coordinate. In this region the flow quantities have expansions of the form: 

u = Jj b (x, T) T eu\M -f i 7 u 2 + ... + c 5 U 5 4" - . * > 
v = t 7 v 2 + . . . + c 5 t? 5 + . . . , 

P = PB + *PlM + t 2 (P2 + P2M) + ... + € 5 p 5 + * , (2.1) 

p — Pb{ x i Y) + CPlM + + . . « + +f 5 ps + . . . , 

T = Tb{ x i Y ) + + . . * + + • - • * 

The terms with suffices M denote mean flow terms which arise because of the properties of 
the nonlinear critical layer. The fundamental travelling wave disturbances have amplitudes 
0(e 7 ) and we write 

with C = cos A r , S = sin X, and where 

d _ 6 d d d _ 6 d 
dx f a °dA + fo’ a ° C °dX' 

The wavenumber a o and the phase speed c<> are taken to be real. After substitution into 
the Navier-Stokes equations the equation governing the stability of the flow is given by, 

L r V 2 = 0 (2.2) 
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where 


d d 

Lr = dY^dY^ ~ ^ + a2 ° PB ^ 


with 


, PB a= d 9U B 1 

* 1 -M 2 ' q dY^ dY (U B - c 0 )’ 


M 


MIoPb(I b - c 0 ) 2 , 


is the compressible inviscid linear stability operator. The problem (2.2) has to be solved 
subject to the boundary conditions, 

V 2 (0) = 0, and V 2 (oo) < oo, 

for no slip at the wall and boundedness in the freestream, together with suitable jump 
conditions across the critical point Y — F c where Ub{Yc) — c o- This point is a regular 
singular point of the differential equation and a Frobenius expansion locally shows that 


V 2 (Y) = D 0 Y,9n(Y ~ Y c ) n + D } 
0 

where 


Y cur - y;r + mY s»( y - y c)") ln (y - y c ) 


go = 0, g\ — 1, Go — 1, G\ — 0, 


A = 


Mpb*&) 

eu B 


dY 


at Y = Y c , 


PB 


dY 


and D 0 , £>i are arbitrary constants. To complete the specification of the problem the jump 
conditions across Y = Y c are necessary. Following the discussion in the introduction the 
appropriate jump conditions for the strongly nonlinear critical layer are that, 


ln(y - Y e ) -» In \Y - Y e \ 


for 


Y <Y C 


The amplitude equation which determines the possible neutral states can be derived as 
in BSG (1983) with suitable modifications to allow for compressibility. This gives the 
equation: 


D Ik \\(!h I !) - ^ ; g'(0) 2 Co / pw g (2 3i 

R oPc(bx) ( Rq + ^ )C ao{Tw - Mlcl)2 { ’ 


where 


. dU b , ol d 2 U B[ D rp „ (v n \ 1 

b i = ~0y\Y c , 2b 2 = —\ Y ', Ri- dY W c , Tw-T B (Y- Q)- pw -. 
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Ro = PB In, C' (1) = iC (,) (l - 2 J) + 2C W J - 2tt7;, 

giji 

( R\ b\ + 2^2 Ro) 

and /x c , i-i w i are the coefficients of viscosity at the critical point and the wall respectively. 
Also the normalization V 2 = Ag(Y) with g{Y c ) — 1 has been used, .4 being a scaled am- 
plitude, and C (1) , G'( 2 \ are the constants appearing in GC equation (4.8). The detailed 

analysis leading to (2.2), (2. 3), which requires consideration of the higher order terms in 
(2.1) and details of the nonlinear critical layer, and the wall layer, is given in Gajjar &: 
Cole (1989b). In summary the eigenvalue problem (2.2) determines Qo as a function of 
the phase-speed Co and the Mach number M <*,, and the corresponding neutral amplitude 
is obtained from (2.3). The scaled amplitude used in BSG differs from that used here but 
is in fact equal to /^(O) and related to A by, 


PA o) = 


cpAff^lO) 

«o (Tw - A/^Co) 


(2.4) 


For (2.3) to hold we require in addition that, 

(^ + ^)C <1)>0 . (2.5) 

Kq 0] 

A special case of the above is when the generalized inflexion point criterion holds, 
namely A = 0. This gives the same wavenumber 00 — a, and phase-speed Co = c, as in 
the linear case, but with the nonlinear neutral amplitude from (2.3). 

For a critical layer to exist in the main part of the boundary layer we require 0 < Co < 
1. Also we restrict ourselves to modes which are subsonic in the freestream (and then V 2 
decays exponentially for large K), so that, 1 — 1 /Mqo < Co <1. 

The basic flow for which most of the results were obtained is the flat plate boundary 
layer flow given by: 


Va = fM, f" + //" = »> /( 0 ) = /*( 0 ) = 0 , /'(<*>) = 1 , 

Tb = 1 + — ~ — Mlo( 1 - / ,2 L f° r insulated walls, 

Tb = 1 + ^ - — M^Ub ~ U%) ~ (1 ~ 7V)(1 - Ub ), with heat transfer. 
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Here q is the Dorodnitsyn-Howarth variable with, 

(2Cx 0 )*ri - f p B dY. 

Jo 

and C is Chapman’s constant. 

The solution(s) of (2.2) were obtained numerically using a shooting technique supple- 
mented by series expansions near Y = Y c . It was found convenient to work in terms of the 
variable 77 . Various checks were made including those for M a 0 — 0, and for values of M ^ 
and c 0 for which the generalized inflexion point criterion was satisfied. 
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3. Results and Discussion 

Some preliminary results obtained are shown in Figures (3-8) It is apparent that there 
are different behaviours depending on whether M oo < 1, M > 1 and whether the flow 
is relatively subsonic or supersonic in some parts of the flowfield. For convenience each of 
these cases is dealt with separately below for insulaled walls and for when there is heat 
transfer. 

(a) Insulated walls. If the flow is everywhere relatively subsonic then M 2 < 1 and this 
gives 1-1 1 Moo < c 0 < c a = + ^fr), whereas if M 2 > 1 then 

Cb = max(c a ,l — 1/Moo) < Co < 1, are the possible values of the phase-speeds. The curves 
depicting these are shown in Fig. 2 where the curve for c = c, is also included. Further 
restrictions on Co arise from ( 2 . 5 ). For c<> > c d, ( c d is the value of c o when C ^ = 0, ) (2.5) 
is satisfied. 

(i) Moo < 1. For this case there is only one neutral mode, the properties of which are 
similar to that of incompressible mode computed in BSG. (Incidentally we find that the 
wavenumbers shown in BSG for Blasius flow should be multiplied by a factor \/2 ). Results 
for Mqo = 0, . 5 , and 1 are shown in Fig. 3. The lower limit for the phase-speed is c 0 = cj 
from (2.5), and as c 0 — * Cd+, the amplitude goes to zero and from (2.3), 


A 


<F + 

Hq Oi 


2 

3 


If Moo is small also this then matches with the corresponding limit in Cole (1989). As 
Co — » 1— the neutral w r avenumber and amplitude increase indefinitely, for fixed values of 
Mqo • This limiting structure is similar to that in Smith et al. (1989) for the corresponding 
incompressible mode, and we find that 


(Tw - Ml 


2 \* - 1 
i ) 3 ^ 7 


where is the asymptotic location of the critical point given by, 


c -j(<r~Th ) 7 

332 
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and 7/1 = 1.2167. 

(ii) M 2 < 1 , Moo > 1* Since the flow is everywhere subsonic there is only one neutral 
mode, but with phase-speeds lying in the restricted range max(c<f,l - 1 /Moo) < c 0 < c a . 
Figure 4 shows the variation of the neutral wavenumber and the neutral amplitude, against 
phase-speed for Afoo — 2.2. As Co approaches 1 — l/M^, the wavenumber goes to zero. We 
find that 

a 0 ~(1 - Ml) 1 ', A ~ (1 - Ml)$ , P 2 (0)^{l-MlrK 

consistent with the numerical results. For Co — * c a — the wavenumber and amplitude 
approach constant values the reason being that g'(0) approaches zero, 

*'( 0 ) - (Tw ~ Mlcl). 

(iii) M 2 > 1. The flow becomes relatively supersonic near the wall and the differential 
equation changes type. As in the linear case, Mack (1984), there exist infinitely many 
neutral modes. The first four were computed and are shown in Figure 5, together with the 
corresponding neutral amplitudes. It can be seen that for certain values of Co, for fixed 
values of M 0 c , the modes after initially decreasing start to increase, and then switch to take 
on the characteristics of the next higher mode. This is also seen in Figures 6(a-c), which 
show the eigenfunctions for different values of c<>. The sudden change in character of one 
mode to the next higher mode can be explained as being due to the presence/superposition 
of two distinct modes, one say, the WKB mode, and the other the ’subsonic’ mode men- 
tioned in (i) above, as can be seen clearly in Fig. 6. An analysis incorporating these ideas 
gives the wavenumber implicitly as, 

a, - a> = «t + \)* + 2Uo- | | F ' (0 ^ ( 1 ° ll)/r(r.), (k =0,1,...), (3.2) 

4 a*<£'(0) 

where 

Ml(Y.) = 1, F(Y ) = f \\-M'\iiy. 

Jo 

The first term in the bracket on the right in (3.2) is the WKB mode, and the second 
describes the additional contribution from the ’subsonic’ mode. The contribution from the 
WKB mode alone in (3.2) is superimposed on the computed values in Figure 5a and the 
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agreement is quite good. The change in character of the mode can be explained by the 
presence of the second term in which ( c is as in (3.1) the location of the critical point. 
The quantity (p comes from solving the reduced eigenvalue problem for the subsonic mode 
but with slightly changes boundary conditions. For (, c J a k << 1, the second term does not 
contribute. As ftfc decreases then at some critical values of c 0 we have ( c /°* - 0(1), and 
then the second term gives a nonzero contribution and a* increases. For ( c / a k >> 1 the 
right side of (3.2) is exactly the same as that for the (khl)th VVKB mode. 

The limiting form (3.2) can be used to give the appropriate behaviour of A, but this is 
not done here. We note that a special case of the above, corresponding to the generalized 
inflexion point, is when Co = c, and «o = a,. Results for this case, for linear theory, similar 
to those described above have been obtained independently also by Cowley & Hall (1989) 
( they actually considered the limit for M «> >> 1 ), and Brown & Smith (unpublished but 
mentioned by Cowley & Hall ). Our results agree with those of Cowley & Hall (1989) for 
Moo >> 1 since from (2.4) a* = (^(l/M^). The linear results for the generalised inflexion 
point mode given in Mack (1984) are different from those obtained here. It is not clear 
what viscosity law Mack has used, and this may be one reason for the difference. We also 
note that we could not obtain the first peak in the plot of the wavenumber against Mach 
number as shown in Fig. (9.2) in Mack (1984) although the second peak near M = 4 was 
obtained. We cannot explain this discrepancy, and in fact our results are entirely consistent 
with the analysis and asymptotes in GC and here, and also Cowley &: Hall (1989). 

Finally for large M <*, the neutral modes are confined to a narrow range of phase-speeds 
close to the freestream value, as can be seen also from Fig. 2. 

(B) Heat transfer Only a few calculations were carried for cases when there is 
heat transfer. Some typical results are shown in Figures 7,8. For fixed values of the wall 
temperature Tw , the results obtained are broadly similar to those described above for the 
insulated case. For small values of c 0 and Moo < 1, the asymptotes given in GC apply, and 
the numerical results agreed well with these. For M 2 > 1 the asymptote (3.2) still holds 
and a comparison with the computed values is shown in Figure 8a, for the WKB mode. 

In summary, we have considered the compressible version of the strongly nonlinear 
critical layer of BBD. Some new neutral modes have been computed and these have different 
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properties depending on whether the flow is relatively subsonic or supersonic in some 
regions. In the latter case there exist infinitely many neutral states. Which if any of 
these can be observed in practical situations or in full numerical computations, is an open 
question. The work could be extended further to consider the next stage when the various 
limits are approached or when >> 1, but in our opinion it would be more useful to 
follow the spatial/ temporal evolution of a linear unstable mode into the nonlinear stages 
of evolution. This would then bring into play some of the ideas discussed here and also in 
GC, on nonlinear critical layers, and address the more important issue of whether any of 
these neutral states are significant or not. Work on this is currently in progress. 
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FIGURE 1, - NONLINEAR CRITICAL LAYER DISTURBANCE STRUCTURE. 
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(a) NEUTRAL WAVENUMBER a Q = a 0 (2x 0 C) AGAINST PHASE-SPEED 
FIGURE 3. - BLASIUS FLOW WITH IN! 


FIGURE 2. - PHASE-SPEEDS FOR BLASIUS FLOW WITH ADIABATIC 
WALL CONDITIONS: C Q = C s IS PHASE-SPEED FOR GENERALIZED 

INFLEXION POINT; c Q = c a IS PHASE-SPEED FOR WHICH M 2 = 1 
AT THE WALL; c 0 = c d IS CURVE FOR C (1) * 0. THE UNLABELLED 
CURVE IS C 0 = 1 - 1/Woo. 
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(b) CORRESPONDING NONLINEAR NEUTRAL AMPLITUDES. 
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(b) CORRESPONDING NONLINEAR 
NEUTRAL AMPLITUDES. 


FIGURE 4. - BLASIUS FLOW WITH INSULATED WALL CONDITION FOR M^ = 2.2. 



(a) NEUTRAL WAVENUMBER AGAINST PHASE-SPEED. (b) CORRESPONDING NONLINEAR NEUTRAL AMPLITUDES FOR FIRST 

FOUR MODES. 

FIGURE 5. - BLASIUS FLOW WITH INSULATED WALL CONDITIONS FOR M^ * 5. 
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(a) NEUTRAL WAVENUMBER AGAINST PHASE-SPEED. (D> CORRESPOND In 

FIGURE 7. - BLASIUS FLOW WITH HEAT TRANSFER FOR M = 0.5 


(b) CORRESPONDING NEUTRAL AMPLITUDE. 


NUMERICAL CALCULATIONS 
WKB ASYMPTOTES (EQ. (3.2)) 
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(a) NEUTRAL WAVENUMBER AGAINST PHASE-SPEED. 
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(b) CORRESPONDING NONLINEAR NEUTRAL AMPLITUDES FOR FIRST 
FOUR MODES. 


FIGURE 8. - BLASIUS FLOW WITH HEAT TRANSFER FOR = 5 AND T w = 3. 
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